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ACQUIRING A FREQUENCY AND PHASE OFFSET 
ESTIMATES USING FREQUENCY DOMAIN ANALYSIS 

GOVERNMENT FUNDING 

The U.S. Government may have certain rights in this 
invention as provided for by the terms of Grant Nos . MUOS 
G8238 TCCPM awarded by the U.S. Navy. 

TECHNICAL FIELD 

This invention relates generally to the field of 
signal processing and more specifically to acquiring a 
frequency and phase offset estimates using frequency 
domain analysis. 
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BACKGROUND 

Certain channel coding techniques allow wireless 
receivers to operate at lower signal-to-noise ratios than 
before. As a result, estimating carrier offsets to track 
signals is more challenging. Some known techniques use a 
preamble or amble to remove data modulation, then 
estimate frequency offsets using phase differentiation 
with respect to an elapsed time period. These 
techniques, however, typically increase noise, which 
reduces performance at low signal-to-noise ratios. Other 
known techniques rely on the use of lengthy continuous - 
wave (CW) or alternating sequence (AS) preambles. These 
techniques, however, typically result in spectral 
interference and are well known for their false 
acquisition. 
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SUMMARY OF THE DISCLOSURE 

According to one embodiment of the present 
invention, determining a frequency and phase offsets 
estimate includes receiving a signal at an offset 
5 estimator- The received signal is zero-padded in the 

time domain to yield a zero-padded signal. A Fourier 
transform of the zero-padded signal is taken to yield a 
transformed signal. The maximum power of the transformed 
signal is established. A frequency offset estimate is 
10 generated based on the maximum power of the transformed 

signal . The frequency offset may be used to estimate a 
phase offset . 

Certain embodiments of the invention may provide one 
or more technical advantages. A technical advantage of 
15 one embodiment may be that a frequency offset is 

estimated using frequency domain analysis. The embodiment 
does not require phase differentiation in the time domain 
analysis, which may increase noise. Another technical 
advantage of one embodiment may be that a short preamble 

2 0 may be used. The embodiment does not require lengthy 

continues wave (CW) or alternating sequence (AS) 
preambles, which typically contribute to undesired 
spectral interference and result in false acquisition. 
Yet another technical advantage of one embodiment may be 
25 that a digital feed- forward (open loop) topology may be 

used. The topology may avoid problems associated with 
feedback methods that deteriorate when used in multipath 
wireless channels. 

Certain embodiments of the invention may include 

3 0 none, some, or all of the above technical advantages. 

One or more other technical advantages may be readily 



ATTORNEY'S DOCKET NO. 
064750 . 0475 



PATENT APPLICATION 



apparent to one skilled in the art from the figures, 
descriptions, and claims included herein. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present 
invention and its features and advantages, reference is 
now made to the following description, taken in 
conjunction with the accompanying drawings, in which: 

FIGURE 1 is a block diagram illustrating one 
embodiment of a system for acquiring frequency and phase 
offsets; and 

FIGURE 2 is a flowchart illustrating one embodiment 
of a method for acquiring frequency and phase offsets 
that may be used with the system of FIGURE 1 . 



ATTORNEY ' S DOCKET NO . 
064750 . 0475 



PATENT APPLICATION 



6 

DETAILED DESCRIPTION OF THE DRAWINGS 

Embodiments of the present invention and its 
advantages are best understood by referring to FIGURES 1 
and 2 of the drawings, like numerals being used for like 
5 and corresponding parts of the various drawings. 

FIGURE 1 is a block diagram illustrating one 
embodiment of a system 10 for acquiring frequency and 
carrier phase offsets. System 10 may acquire frequency 
and carrier phase offsets of a modulated signal using 

10 frequency domain analysis applied to a data-aided 

technique. System 10 may track the signal using a 
decision-aided technique. 

According to the illustrated embodiment, system 10 
includes an input 20, one or more pre-processing modules 

15 22, an offset estimator 24, and one or more error 

correction modules 26 coupled as shown. According to one 
embodiment of operation, input 20 receives a signal 
communicating symbols. Pre-processing modules 22 perform 
pre-processing operations on the received signal. Offset 

20 estimator 24 estimates a frequency offset and a carrier 

phase offset of the signal using frequency domain 
analysis. Error correction modules 26 estimate residual 
error and correct the received signal in accordance with 
the frequency offset estimate, the phase offset estimate, 

25 and the residual error estimate. 

According to one embodiment, a signal may be 
received at a receiver antenna system. A radio frequency- 
to- intermediate frequency converter and an analog- to- 
digital converter may convert the signal. The converters 

30 may then feed a sampled digital signal to input 20. 

A signal refers to an electrical quantity, such as a 
current or voltage, that may be used to convey symbols 
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through a channel. A signal may comprise, for example, a 
coded signal. The signal may include packetized 

information. A packet 2 7 of the signal may have any 
suitable format and size. Data 28 of packet 27 may be 
5 preceded by a pseudo-noise (PN) code 29 that may be less 

than approximately 10% of the total packet length- PN 
code 2 9 may be known prior at the receiver, then 
exploited by offset estimator 24 using multiplier 40. 

Pre-processing modules 22 perform pre-processing 

10 operations on the received signal such as removing 

channel gain from the signal, removing sample timing from 
the signal, decimating the signal, other suitable pre- 
processing operations, or any combination of the 
preceding- Pre-processing modules 22 of this embodiment 

15 include an automatic gain controller (AGC) 30, a delay 

32, a multiplier 34, and a decimator 36 coupled as shown. 

Automatic gain controller 32 estimates the channel 
gain of the received signal. Automatic gain controller 
3 0 may comprise one or more detectors operable to 

2 0 maintain a constant amplitude of a signal by reducing or 

increasing the gain of the signal in accordance with the 
strength of the signal- Delay 32 delays the received 
signal for a suitable amount of time such that the signal 
and the channel gain estimate for the signal arrive at 
25 multiplier 34 at substantially the same time- Delay 32 

may comprise a delay circuit, such as a first-in-first- 
out (FIFO) circuit that introduces a time delay- 
Multiplier 34 removes the channel gain from the received 
signal- Decimator 3 6 reduces the rate of the received 

3 0 signal. 

Offset estimator 24 estimates a frequency offset, a 
carrier phase offset, or both a frequency offset and a 
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carrier phase offset of the received signal using 
frequency domain analysis. A frequency offset refers to 
the difference between a transmitter frequency and a 
reference frequency, and a phase offset refers to the 
5 difference between the phase of the received signal and a 

reference phase, where the difference is typically 
introduced by the channel through which the received 
signal has passed. Offset estimator 24 includes a 
multiplier 40, a transformer 42, a frequency offset 

10 estimator 44, and a phase offset estimator 46 coupled as 

shown. Multiplier 40 removes the modulation of the 
received signal using a pseudo random code 48 to yield a 
signal with a residual tone that corresponds to the 
actual frequency offset . 

15 Transformer 42 pads a predetermined window of time 

that is thought to include PN code 2 9 with zeros and then 
performs a Fourier transform on the zero-padded signal to 
yield a frequency transformed signal. Zero-padding refers 
to adding one or more zero-valued samples to a signal . 

2 0 Any suitable number of zero- valued samples may be added 

to any suitable portion of a signal, such as the 
beginning, middle, or both the beginning and middle of 
the signal. A Fourier transform may refer to a discrete 
Fourier transform, such as a fast Fourier transform, used 

25 for frequency analysis of discrete signals. Although 

transformer 42 is shown performing the zero-padding and 
the fast Fourier transform, any number of modules may 
perform these operations. Frequency offset estimator 44 
estimates the frequency offset of the transformed signal, 

30 and provides the frequency offset estimate to phase 

offset estimator 46. Phase offset estimator 46 estimates 
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the phase offset from the same transformed signal and .the 
frequency offset estimate. 

Error correction modules 26 estimate the residual 
error and correct the received signal for the frequency 
5 offset, the phase offset, and the residual error. Error 

correction modules 2 6 include a numerically controlled 
oscillator (NCO) 50, a multiplier 54, a data decoder 56, 
a phase detector 58, a delay 60, and a loop filter 62 
coupled as shown. Error correction modules 2 6 are coupled 

10 to a delay 52 as shown. 

Delay 52 delays the received signal for suitable 
amount of time such that the received signal and 
correction information from numerically controlled 
oscillator 50 arrive at multiplier 54 at substantially 

15 the same time. Delay 52 may comprise a delay circuit, 

such as a FIFO circuit that introduces a time delay. 
Multiplier 54 applies the correction information to the 
received signal to yield a corrected signal. Decoder 56 
decodes the corrected signal to yield a decoded signal 

20 having the symbols of the received signal. Delay 60 

delays the corrected signal by a single symbol period T 
such that the delayed signal and the decoded signal 
arrive at phase detector 58 at substantially the same 
time. Delay 60 may comprise a delay circuit, such as a 

25 FIFO circuit that introduces a time delay . 

Phase detector 58, along with loop filter 62 and 
numerically controlled oscillator 50, form a phase-locked 
loop. Phase detector 58 compares the received delayed 
signal corrected using estimates from offset estimator 24 

30 with the decoded symbols estimate residual error. Phase 

detector 58 may comprise a detector for phase modulation 
and frequency modulation, and may generate an estimate 
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corresponding to the difference in phase between the 
signals. Loop filter 62 determines an error correction 
from the residual error estimate. 

Numerically controlled oscillator 50 receives 

5 frequency offset estimates , determines corrections in 

accordance with the estimates, and sends correction 
information to multiplier 54. For example, numerically 
controlled oscillator 50 receives the phase offset 
estimate from phase detector 58, determines phase offset 

10 corrections 0^ from estimator 24 and residual errors 

from the loop filter 62, and sends correction information 
to multiplier 54. Numerically controlled oscillator 50 
may include one or more circuits operable to generate a 
digital oscillating signal. As described previously, 

15 multiplier 54 applies the correction information to the 

received signal . 

System 10 may be used for wireless receivers that 
operate at a low signal-to-noise ratio. These wireless 
receivers may use any suitable channel coding technique, 

20 such as Reed- Solomon coding or turbo coding for constant 

envelope bandwidth efficient modulations such as 
continuous phase frequency shift keying (CPFSK) . 

Alterations or permutations such as modifications, 
additions, or omissions may be made to system 10 without 

25 departing from the scope of the invention. System 10 may 

have more, fewer, or other modules. Moreover, the 
operations of system 10 may be performed by more, fewer, 
or other modules. For example, the operations of 

transformer 42 may be performed by more than one module. 

30 Additionally, operations of system 10 may be performed 

using any suitable logic comprising software, hardware, 
other logic, or any suitable combination of the 
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preceding. As used in this document, "each" refers to 
each member of a set or each member of a subset of a set . 

FIGURE 2 is a flowchart illustrating one embodiment 
of a method for acquiring frequency and carrier phase 
5 offsets that may be used with the system of FIGURE 1. 

Although the method is described as used with system 10 
of FIGURE 1, the method may be performed using any 
suitable system for acquiring frequency and carrier phase 
offsets . 

10 The method begins at step 100, where a signal is 

received. According to the illustrated embodiment, input 
20 receives a signal that has traveled through a channel. 
A packet of the signal may have, for example, a preamble 
of L=50 symbols and a total coded packet size of Ltot = 

15 1024 symbols with Minimum Shift Keying (MSK) modulations. 

The preamble symbols {c^} may be less than approximately 

five or ten percent of the total packet size. 

The signal is pre-processed at step 102. For 
example, the signal may be sampled by a receiver analog- 
20 to-digital (A/D) converter sampling at a rate of l/T^ , 

which is an integer multiple of the symbol rate l/T . The 
sampled signal y^^ = y{kT) may be given by Equation (1) : 



25 



yk - Sk^k^ ^'^k (1) 



where k represents the sample Ts index, represents the 

transmitted complex data symbol sample. A/* represents the 
frequency offset, gj^ represents the channel gain 
magnitude sample, 0^^ represents the phase shift 
3 0 introduced by the channel, x represents the sample timing 
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error, and n^^ represents the complex additive white 

Gaussian noise (AWGN) . 

According to the illustrated embodiment, automatic 
gain controller 3 0 estimates the channel gain g^^ . 
5 Multiplier 34 removes channel gain g^^ from signal y{kT^) 

delayed by delay 32. The sample timing may also be 

removed. To reduce computational complexity, decimator 
36 may decimate signal y{kT^) to a single sample per 
symbol to yield signal y{kT) . 

10 The transmitted modulations of the signal are 

removed at step 104, The transmitted modulations for 
symbol sample may be removed by multiplying signal y^^ 
by a locally generated conjugate replica 48 of preamble 
symbols {cD , where * indicates a complex conjugate. 

15 According to the illustrated embodiment, multiplier 40 

receives pseudo- random code a^^* 4 8 and signal to yield 
unmodulated signal . If the sample timing {T^-x) is 
accurate, signal may be given by Equation (2) : 

20 z,=3',c;=cXe>(^'^-«'> + «,c: (2) 

Since clUf^-l , the baseband signal given by Equation (2) 
may also be described by Equation (2 • ) : 

2 5 z, =e^'('"^^"®*) + «,c; (2«) 

The signal is filtered at step 108 to pass through 
the maximum expected frequency offset. For example, a 
linear bandpass filter (BPF) may have a bandwidth that is 
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only wide enough to pass through the maximum expected 
frequency offset BW^pj, = Af^^ . The filtering may reduce 
the probability of false frequency offset detection due 
to noise or spurious receiver emissions. 
5 The signal is zero-padded at step 112. According to 

the illustrated embodiment, transformer 42 zero-pads 
signal . Any suitable number of zero-valued samples 
may be added to the signal to zero-pad the signal in any 
suitable manner. According to one embodiment, to 

10 generate a fast Fourier transform of N samples for L 

symbols, N-L zero-valued complex samples may be appended. 
For example, to generate a fast Fourier transform of 
N=1024 samples for L=50 symbols, N-Lj=974 zero-valued 
complex samples may be appended. 

15 A fast Fourier transform (FFT) of the signal is 

taken at step 116. According to the illustrated 

embodiment, transformer 42 takes a fast Fourier transform 
of signal . The resulting transformed signal may be 
described by Equation (3) : 



20 



Z(K) = Z(^'^''^'^'*^ (3) 



*=0 



where k represents the frequency bin index, k = fjN 
represents the FFT bin resolution, and represents the 
25 sampling frequency, which in this case is equal to the 

symbol rate (l/T) . Any suitable number N of FFT bins may 
be used, for example, N=1024, where k ranges from k = 0 to 
k = N-l, The N bins may include a single bin with maximum 

enerqy A: = k , which corresponds to the frequency offset 
max ^ ^ 

30 estimate . Because a finite number of bins are used. 
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the actual frequency offset 6f may not coincide with the 
a bin, but may spread over multiple bins including 
A: = k 

max 

The frequency offset is estimated at step 120. 
According to the illustrated embodiment, frequency offset 
estimator 44 determines frequency offset estimate A/* by 
locating the FFT bin corresponding to the maximum power. 
Frequency offset estimate may be given by Equation 

(4) : 



A/ = max|Z(K)| (4) 



Using the estimate obtained by Equation (4) , the 
carrier phase offset is estimated at step 124. According 
15 to the illustrated embodiment, phase offset estimator 4 6 

determines carrier phase offset estimate 9^ according to 

Equation (5) : 



=arg(Z(K_))-27cA/r (5) 

20 

where arg (x) represents the argument of complex value x 
computed using the arctangent function, K^ax is the index 
of the FFT bin corresponding to the maximum power, and k 
in this case refers back to the symbol index. The 
2 5 disturbance terms in Equations (4) and (5) may be 

determined by expanding Equation (3) as Equation (3 • ) : 



ZiK)=Ze^ + " (3-) 



disturbance 
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The offset estimates based on Equations (4) and (5) may 
have fewer disturbance terms than estimates obtained 
using conventional data-aided methods. Moreover, the 
zero-padding performed at step 112 may improve offset 
estimates by averaging the estimates over a longer 
interval N than the interval used in conventional data- 
aided methods . 

The signal is corrected at step 128. According to 
the illustrated embodiment, numerically controlled 
oscillator (NCO) 50 receives offset estimates from offset 
estimator 24 and a residual error correction from loop 
filter 62. Numerically controlled oscillator 50 sends 
correction information to multiplier 54 to correct the 
offsets and residual error of signal delayed by delay 

52 to yield signal y^^ given by Equation (6) : 

j;,=a,e^^^*^^^^^"^^>^+«, (6) 

20 where rijc is the resulting additive white Gaussian noise 

AWGN, k represents the symbol index, 5t represents any 
timing residual error, 5/ represents the frequency offset 
residual error, and 59 represents the carrier phase 
offset residual error. The frequency offset residual 

2 5 error and the carrier phase offset residual error result 

from inaccuracies of offset estimators 44 and 46. The 
signal is decoded at step 130. According to the 
illustrated embodiment, decoder 56 decodes the signal to 
yield symbols . 

3 0 The residual error is estimated from the decoded 

signal at step 132 . The carrier recovery parameters 



10 
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based on Equations (4) and (5) may have some residual 
errors due to FFT bin resolution. The undesired phase 
rotation error ^n^k+^'^^iT^^)) removed using decision- 

aided tracking methods. If the network timing is 

5 accurate, the timing residual error 5x may be assumed to 

be negligible, and so only the error offset term e^^^''*^'^^^ 
is tracked and removed. According to the illustrated 
embodiment, phase detector 58 obtains an error e^ik) by 

multiplying the conjugated symbol from decoder 54 and 

10 its corresponding input sample Uj^e^^^"'^^''^^^^ + stored in 

delay 60, and then taking the imaginary part of the 
result, yielding error signal e^ik) given by Equation (7): 

e,(A:) = Im{a;a,e^<-*-->^^>+n,a;} (7) 



15 



20 



Assuming that the transmitted symbols are designed such 
that |<2jj||<2jt| = 1 / and ignoring the noise contribution by the 

term n^al , the error signal e^ik) may approximated using 
small angle approximation by Equation (8) : 

e^{k) = de, ^InbfkT (8) 



where k is the symbol index. 

Loop filter 62 and numerically controlled oscillator 
25 50 operate on error signal ^(^) to form a phase locked 

loop (PLL) that is used to track the final carrier 
offsets. Loop filter 62 calculates a signal given by 
Equation (9) : 
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(9) 



where the constants and represent the proportional 

and integrator loop coefficients, respectively. The 
5 error signal described by Equation (9) is then sent to 

numerically controlled oscillator 50, which uses the 
error signal to correct the signal using multiplier 54. 

Depending on the desired tracking loop dynamics, the 
maximum frequency residual error 5/^ in Equation (7) that 
10 the loop can cope with is bounded by 5/^^<5^, where is 

the tracking loop bandwidth that is largely dependent on 
and K. in Equation (9) . The value of 8/^^ may be used 

to determine the value of N of Equation (4) . 

If decoding is to continue at step 136, the method 

15 returns to step 128 to correct the signal. If decoding 

is not to continue at step 136, the method terminates. 

According to one example, for a packet size of 1,024 
Minimum- Shi ft Keying symbols, the method may achieve a 
perfect carrier offset estimation at EjN^=—2dB , where 

20 Es/No represents the ratio of energy per symbol over the 

noise within one symbol duration, while using only 50 
Minimum-Shift Keying PN symbols. Such performance may be 
suitable for receivers that are required to receive turbo 
coded satellite signals or land based long-range line-of- 

25 sight communication networks requiring at least 

EjN^=2dB for bit-error rate of 10'^. 

Alterations or permutations such as modifications, 
additions, or omissions may be made to the method without 
departing from the scope of the invention. The method 

30 may include more, fewer, or other steps. Additionally, 
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steps may be performed in any suitable order without 
departing from the scope of the invention. 

Certain embodiments of the invention may provide one 
or more technical advantages . A technical advantage of 
one embodiment may be that a frequency offset is 
estimated using frequency domain analysis. The embodiment 
does not require phase differentiation time domain 
analysis, which may increase noise. Another technical 
advantage of one embodiment may be that a short preamble 
may be used. The embodiment does not require lengthy 
continues wave (CW) or alternating sequence (AS) 
preambles, which typically contribute to undesired 
spectral interference and result in false acquisition. 
Yet another technical advantage of one embodiment may be 
that a digital feed- forward (open loop) topology may be 
used. The topology may avoid problems associated with 
feedback methods that deteriorate when used in multipath 
wireless channels . 

While this disclosure has been described in terms of 
certain embodiments and generally associated methods, 
alterations and permutations of the embodiments and 
methods will be apparent to those skilled in the art. 
Accordingly, the above description of example embodiments 
does not define or constrain this disclosure. Other 
changes, substitutions, and alterations are also possible 
without departing from the spirit and scope of this 
disclosure, as defined by the following claims. 



